Demonstration of three-photon de Broglie wavelength by projection measurement 
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Two schemes of projection measurement are realized experimentally to demonstrate the de Broglie 
wavelength of three photons without the need for a maximally entangled three-photon state (the 
NOON state). The first scheme is based on the proposal by Wang and Kobayashi (Phys. Rev. A 71, 
021802) that utilizes a couple of asymmetric beam splitters while the second one applies the general 
method of NOON state projection measurement to three-photon case. Quantum interference of 
three photons is responsible for projecting out the unwanted states, leaving only the NOON state 
contribution in these schemes of projection measurement. 
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I. INTRODUCTION 

Photonic de Broglie wavelength of a multi-photon state 
is the equivalent wavelength of the whole system when all 
the photons in the system act as one entity. Early work 
by Jacobson et al utilized a special beam splitter that 
sends a whole incident coherent state to either one of the 
outputs thus creating a Schodinger-cat like state. The 
equivalent de Broglie wavelength in this case was shown 
to be A/ (n) with (n) as the average photon number of the 
coherent state. Such a scheme can be used in precision 
phase measurement to achieve the so called Heisenberg 
limit 0jSSIB] 0i V( n ) m P nase uncertainty. 

Perhaps, the easiest way to demonstrate the de Broglie 
wavelength is to use the maximally entangled photon 
number state or the so called NOON state of the form 

ma 

|jvooivr> = -L(|jv> 1 |q> 2 + |o>i|Jv> a ), (i) 

where 1,2 denote two different modes of an optical field. 
The N photons in this state stick together either all in 
mode 1 or in mode 2. Indeed, if we recombine modes 1 
and 2 and make an N-photon coincidence measurement, 
the coincidence rate is proportional to 

R N oc 1 + cos(27riVA/A), (2) 

where A is the path difference between the two modes 
and A is the single-photon wavelength. Eq.© shows an 
equivalent de Broglie wavelength of X/N for N photons. 

NOON state of the form in Eq.JlJ for TV = 2 case 
was realized with two photons from parametric down- 
conversion, which led to the demonstrations of two- 
photon de Broglie wavelength yjQ- For N > 2, however, 
it is not easy to generate the NOON state. The diffi- 
culty lies in the cancellation of all the unwanted terms 
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of \k,N — k) with k ^ 0, N in an arbitrary N-photon 
entangled state of 

N 

\$n) =^c k \k,N-k). (3) 

A number of schemes have been proposed |jj El O] and 
demonstrated 0, which were based on some sort of 
multi-photon interference scheme for the cancellation. 

Without exceptions, the proposed and demonstrated 
schemes 0, E3 El El EI for the NOON state genera- 
tion rely on multi-photon coincidence measurement for 
revealing the phase dependent relation in Eq.(0l. Since 
coincidence measurement is a projective measurement, 
it may not respond to all the terms in Eq.Q. Indeed, 
Wang and Kobayashi E3 applied this idea to a three- 
photon state and found that only the NOON state part 
of Eq. J2I contribute to a specially designed coincidence 
measurement with asymmetric beam splitters. The coin- 
cidence rate shows the signature dependence in the form 
of Eq.((5J) on the path difference for the three- photon de 
Broglie wavelength. Another projective scheme was re- 
cently proposed by Sun et al [l5| and realized experimen- 
tally by Resch et al 16] for six photons and by Sun et al 
|l7j for four photons. This scheme directly projects an 
arbitrary N-photon state of the form in Eq.|j3J) onto an 
N-photon NOON state and thus can be scaled up to an 
arbitrary N-photon case. 

In this paper, we will apply the two projective schemes 
to the three-photon case. The three-photon state is pro- 
duced from two pairs of photons in parametric down- 
conversion b y g ating on the detection of one photon 
among them |l8j . We find that because of the asymmet- 
ric beam splitters, the scheme by Wang and Kabayashi 
E3 has some residual single-photon effect under less per- 
fect situation while the NOON state projection scheme 
cancels all lower order effects regardless the situation. 

The paper is organized as follows: in Sect JL we will 
discuss the scheme by Wang and Kobayashi [1J| and its 
experimental realization. In Sect. Ill, we will investigate 
the NOON state projection scheme for three-photon case 
and implement it experimentally. In both sections, we 
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will deal with a more realistic multi-mode model to cover 
the imperfect situations. We conclude with a discussion. 



II. PROJECTION BY ASYMMETRIC BEAM 
SPLITTERS 

This scheme for three-photon case was first proposed 
by Wang and Kobayashi to use asymmetric beam 
splitter to cancel the unwanted |2, 1) or |1,2) term and 
is ageneralization of the Hong-Ou-Mandel interferometer 
[7|,ll,Ei f° r two-photon case. But different from the two- 
photon case, the state for phase sensing is not a three- 
photon NOON state since only one unwanted term can 
be cancelled and there is still another one left there. So a 



special arrangement has to be made in the second beam 
splitter to cancel the contribution from the other term. 
The following is the detail of the scheme. 



A. Principle of Experiment 

We first start with a single mode argument by Wang 
and Kobayashi. The input state is a three-photon state 
of |2) a |l)b. The three photons are incident on an asym- 
metric beam splitter (BS1) with T ^ R from two sides 
as shown in Fig.l. The output state can be easily found 
from the quantum theory of a beam splitter as [2(1 f21| 



|BSl) ouf = V3T 2 i? |3 C , d ) + V3TR 2 |0 C , 3 d ) + VT{T - 2R)\2 C: l d ) + Vi?(P - 2T)\1 C , 2 d ). (4) 

I 



When R = 2T = 2/3, the |l c ,2d) term disappears from 
Eq.© due to three-photon interference and Eq.© be- 
comes 

|BSl) out = ^|3 c ,0 d ) + ||0 G ,3 d ) - ^|2 C) Li). (5) 

But unlike the two-photon case, the |2 C , Id) term is still 
in Eq.© so that the output state is not a NOON state 
of the form in Eq. © . 

Now we can arrange a projection measurement to take 
out the |2 C , Id) term in Eq.©. Let us combine A and B 
with another beam splitter (BS2 in Fig.l) that has same 
transmissivity and reflectivity (R = 2T = 2/3) as the 
first BS (BS1). According to Eq.©, |2 C , l d ) will not con- 
tribute to the probability P 3 (l e , 2/). So only |3 C , d ) and 
|0 C , 3d) in Eq.© will contribute. The projection mea- 
surement of P3(l e ,2^-) will cancel the unwanted middle 
terms like |2 C , Id) from Eq.©. Although the coefficients 
of |3 c ,0d) and |0 c ,3d) in Eq.© are not equal, their con- 
tributions to i-3(l e ,2/) are the same after considering 
the unequal T and R in BS2. So the projection measure- 
ment of i"3(l e , 2/) is responsive only to the three-photon 
NOON state. Use of an asymmetric beam splitter for the 
cancellation of |2 e , Id) was discussed by Sanaka et al in 
Fock state filtering |1S| . 

The above argument can be confirmed by calculating 
the three-photon coincidence rate p3(l e ,2/) directly for 
the scheme in Fig.l, which is proportional to [22|: 

P3(le,2/) = (2 a ,l 6 |et/t 2 / 2 e|2 a ,l fe ), (6) 

with 

e={c + e^V2d)/V3, 
/ = (e^-^c)A/3, U 




FIG. 1: Arrangement of asymmetric beam splitters of a three- 
photon interferometer for the demonstration of the three- 
photon de Broglie wavelength. 



where we introduce a phase ip between A and B. But for 
the first BS, we have 

c= (a + V2b)/V3, 

d= (b-V2a)/V3. 1 ' 

Substituting Eq.Q into Eq.© with Eq.©, we obtain 

P 3 (l e ,2/) = (2 a ,l b |et/t 2 / 2 e|2 a ,l b ) 
32 

= — (l + cos3<^), (9) 

which has a dependence on the path difference A = 
(p\/2ir that is same as in Eq.© but with N = 3. 

B. Experiment 

Experimentally, asymmetric beam splitters are real- 
ized via polarization projections as shown in Fig. 2 and 
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FIG. 2: Experimental setup. HWP1 and HWP2 is set for 
different measurement, PS is the phase shifter between H- 
photon and V-photon. Insets: (a) arrangement with asym- 
metric beam splitters by polarization beam splitter (PBS) 
and (b) the three-photon NOON state projection. 



its inset (a), where a three-photon state of |2#, ly) is in- 
cident on a combination of two half wave plates (HWP1, 
HWP2) and a phase retarder (PS). The first half wave 
plate (HWP1) rotates the state |2#, ly) by an angle a 
to \2 a , l t ) with 



an = a cos a + b sin a, 
ay — b cos a — a sin a, 



(10) 



where cos a = VT = 1/V3 is the amplitude transmissiv- 
ity of the asymmetric beam splitter. Eq. l|10|) is equivalent 
to Eq.|J7||. The phase retarder (PS) introduces the phase 
shift <p between the H and V polarization. The second 
half wave plate (HWP2) makes another rotation of the 
same angle a for the two phase-shifted polarizations and 
the polarization beam splitter finishes the projection re- 
quired by Eq.(|HJ|. 

The three-photon polarization state of |2#, ly) is pre- 
pared by using two type-II parametric down conversion 
processes shown in Fi g. 2. This scheme was first con- 
structed by Liu et al |23| to demonstrate controllable 
temporal distinguishability of three photons. When the 
delay between the two H-photons is zero, we have the 
state of \2h, ly)- In this scheme, two /3-Barium Borate 
(BBO) crystals are pumped by two UV pulses from a 
common source of frequency doubled Ti: sapphire laser 
operating at 780 nm. The H- photon from BBOl is cou- 
pled to the H-polarization mode of BB02 while the other 
V-photon is detected by detector D and serves as a trig- 
ger. The H- and V-photons from BB02 are coupled to 
single-mode fibers and then are combined by a polar- 
ization beam splitter (PBS). The combined fields pass 
through an interference filter with 3 nm bandwidth and 
then go to the assembly of HWP1, PS, and HWP2 to 
form a three-photon polarization interferometer. There 
are two schemes of projection measurement. In this sec- 
tion, we deal with the first scheme in inset (a) of Fig. 2, 
which consists of a PBS for projection and three detec- 
tors (A, B, C) for measuring the quantity P3(l e ,2/) in 
Eq.© by three-photon coincidence. To realize the trans- 
formation in Eqs.Q |HJ), HWP1 and HWP2 are set to 
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FIG. 3: Experimental result for projection measurement with 
asymmetric beam splitters. The data is least-square-fitted to 
Pio(l + V 3 cos3v9 + Vi cos if) with V 3 = 85% and Vi = 5% 
after background subtraction. 



rotate the polarization by a — cos' -!(1/V3) = 54.7°. 
In order to obtain an input state of \2n, ly) to the in- 
terferometer, we need to gate the three-photon coinci- 
dence measurement on the detection at detector D. In 
this way, we ensure that the two H-photons come from 
two crystals separately. Otherwise, we will have an in- 
put state of \2h, 2y). The delay (AX#) between the two 
H-photons from BBOl and BB02 as well as the delay 
(ATy) between the H- and V-photons are adjusted to 
insure the three photons are indistinguishable in time. 
This is confirmed by the photon bunching effect of the 
two H-photons |23j and a g eneralized Hong-Ou-Mandel 
effect for three photons |18j . 

Four-photon coincidence count among ABCD detec- 
tors is measured as a function of the phase shift ip. 
The experimental result after background subtraction is 
shown in Fig. 3. The data is gathered in 100 sec for 
each point and the error bars are one standard deviation. 
Backgrounds due to three and more pairs of photons are 
estimated from single and two-photon rates to contribute 
1.2/sec on average to the raw data and are subtracted. 
The data clearly show a cos3(yS dependence except the 
unbalanced minima and maxima, which indicates an ex- 
tra cos(y9 dependence. Indeed, the data fit well to the 
function 



Pi = P40 (1 + V 3 cos 3ip + Vi cos ip) 



(11) 



with P 40 = 184, V 3 = 85% and Vi = 5%. The x 2 of the 
fit is 30 and is comparable to the number of data of 25, 
indicating a mostly statistical cause for the error. 

The appearance of the cos ip term in Eq. (|llfl is an indi- 
cation that the cancellation of the |2 C , 1^) and |l c , 2^) is 
not complete in Eqs.QJIHl) an d the residuals mix with the 
|3 C , 0d) and |0 C , 3d) terms to produce the cos <p term. This 
imperfect cancellation is not a result of the wrong T, R 
values but is due to temporal mode mismatch among the 
three photons in the input state of |2# , ly). To account 
for this mode mismatch, we resort to a multi-mode model 
of the parametric down-conversion process. 
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C. Multi-mode analysis of three-photon 
interferometer with asymmetric beam splitters 

We start by finding the multi-mode description of 
the quantum state from two parametric down-conversion 
processes. Since the first one serves as the input to the 
second one, we need the evolution operator for the pro- 
cess, which was first dealt with by Ghosh et al [24[ and 
later by Ou j2f| and by Grice and Walmsley [2(| . In gen- 
eral, the unitary evolution operator for weakly pumped 
type-II process is given by 

U =1 + 7] / duJiduj2^{uJi,ui2)a ji H (uj 1 )ay(uj2), (12) 



where V\ and V denote the two non-overlapping verti- 
cal polarization mode from the first and second crystals, 
respectively. Here we only keep the four-photon term. 
Although there are other four-photon terms in the \^f 2 ) 
state corresponding to two-pair generation from one crys- 
tal alone, they won't contribute to what we are going to 
calculate. So we omit them in Ea. (|14fl . 

The field operators at the four detectors are given by 

E A (t) = nE H (t) + P iE v (t), _ 
E B (t)=[T 2 E v (t)+p 2 E H (t)]/V2 + ... = Ec(t), [ ' 

with n = (1 - 2e^)/3,r 2 = (e^ - 2)/3,p x = -p 2 = 
•\/2(l + e lip )/3. Here we used the equivalent relations in 



where rj is some parameter that is proportional to the 
pump strength and nonlinear coupling. For simplicity 
without losing generality, we assume the two processes 
are identical and are governed by the evolution opera- 
tor in Eci. (|12fl . Furthermore, we assume the symmetry 
relation $(cji,a->2) = $(aj2,wi) which is in general not 
satisfied but can be achieved with some symmetrizing 
tricks |27],|28|. So for the first process, because the input 
is vacuum, we obtain the output state as 



(13) 



(14) 



Eqs.J7||HJ to establish the connection between the field 
operators Ea, Eb, Ec and Eh,Ev and for detector D, 
we have 

E D (t) = E Vl (t), (16) 

with 

E k (t) = -L / dLoa k {Lo)e- wt . (k = H, V, Vi) (17) 

V 27T J 

The four-photon coincidence rate of ABCD is propor- 
tional to a time integral of the correlation function 



I 



I'J'i) = U\vac) = \vac) + r\ J dujiduj2^{i0i, U2)a H {uJi)a\ / {uj2)\vac) . 
The second crystal has the state of as its input. So after the second crystal, the output state becomes 
|*2> = U\^i) = ■■■ + rj 2 I du 1 du2<^[<^2®(u 1 ,U2)§{oj[,U2)a) H (u[^^^ 



r 



T (A) {t U t2MM) = (*2|Sb( < 4)^c7( < 3)^(i2)^(tl)^(il)^B(t2) J Bc(<3)^(i4)|* 2 ). (18) 

It is easy to first evaluate EA(ti)EB(t2)Ec(t 3 )EE,{ti): 

EA(ti)E B (t 2 )Ec(t3)E D (U) = [(HHV + HVH)DT 1 T2P2 + VHHDp 1 p 2 2 ]/2 + ..., (19) 

where H = Eh, V = Ey, D = Ed for short and we keep the time ordering of txt 2 t 3 . We also drop five terms that 
give zero result when they operate on | - It is now straightforward to calculate the quantity E A (ti)EB(t2)Ec(t 3 ) 
ED{ti)\^2) , which has the form of 

2 

E A {t 1 )E B {t 2 )E c {h)E D {t A )\^ 2 ) = |-{ [G(h,t 2 ,t 3 , U) + G(t 2 , *i, t 3 , h) + G{t x ,t 3 , t 2 ,t 4 ) + G{t 3 , t u t 2 , U)] T\T 2 p 2 

+ [G{t 2 ,t 3 , + G(t 3 , t 2 ,t 1 , U)} Pl p 2 2 }\vac) , (20) 
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where 



G(h,t 2 ,t 3 ,U) = g{ti,t 3 )g(t2,U) with g{t,t') = ^- J dLu 1 dLu 2 <S>{Lo U LU 2 )e-^ lt+ " 2t ' ) . 
Substituting Eq. (|20|) into Eq. l|18fl and carrying out the time integral, we obtain 

y duj 1 duj 2 duj' 1 duj 2 [$(uj 1: u 2 )$(u[,u 2 ) + $(u 1 ,u[)$(uj 2 ,uj 2 )](t 1 t 2 p 2 + p 1 pl) 



(21) 
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+ 2$(uj[ , CJ 2 )$(W1, U 2 )TlT 2 p 2 



(22) 



With Ti, pi,T 2 , p2, we can further reduce Ea. (|22|l to 
2M 4 (17.4+7£) 



P 4 oc 
where 

and 



243 



1 + V3 cos 3<p + Vi cos ip 



(23) 



8Q4 + 2£) 9(4 -£) 



174 +7£ 



174 +7£' 



(24) 



4 = / dwid^dwW^l^l^ii ^2)^(^1, u' 2 )\ 2 (25) 

dwi(iw2rftL | J(ia;2 < i ) *(wi, W2) < I > *(^i, w 2 ) 

x$(wi,«a)$(w ls wi). (26) 

In deducing Eas. (|22ll26*|l . we used the symmetry relation 
$(wi, w 2 ) = ^(^2,^1). 

Obviously, when A = £, Eq. l|23|) completely recov- 
ers to Eq.©. In practice, we always have A > £ by 
Schwartz inequality. When £ < A, Ea. l23fl has the 
same form as Ea. l|ll|) . indicating that the multi-mode 
analysis indeed correctly predicts the imperfect cancella- 
tion of the \2 H ,lv) and \2 V , Iff) terms in Eqs.Q EJ. 
If we use the experimentally measured V3 and Vi in 
Eq. (|24fl . we will obtain two inconsistent values of E/A: 
(£/A) 3 = 0.65 and (£/A)i = 0.87. The discrepancy 
is the result of the break up of the symmetry relation 
of <f>(u>i,uj 2 ) = $(u> 2 ,lui) for type-II parametric down- 
conversion, which is reflected in the less-than-unit visi- 
bility of the two-photon interference. This imperfection 
can be modelled as spatial mode mismatch and approxi- 
mately modifies Ea. (|24|l as 



V, = V. 



MA + 2£) _ 9(4 -£) 
17A + 7E' Vl - Vl 17A + 7£> (2?) 



where V\ is the equivalent reduced visibility in single- 
photon interference due to spatial mode mismatch. With 
the extra parameter v\ in Eq. (|27|l , we obtain a consistent 
(£ /A) = 0.86 with vi = 0.96. 



III. NOON STATE PROJECTION 
MEASUREMENT 

The projection measurement discussed in the previous 
section relies on the cancellation of some specific terms 
and therefore cannot be applied to an arbitrary photon 
number. In the following, we will discuss another pro- 
jection scheme that can cancel all the unwanted terms at 
once and thus can be scaled up. 



A. Principle of experiment 

The NOON-state projection measurement scheme was 
first proposed by Sun et al and realized by Resch et 
al 16] for six-photon case and by Sun et al [l7| for the 
four-photon case. Since it is based on a multi-photon 
interference effect, it was recently used to demonstrate 
the temporal distinguishability of an N-photon state p9l 
l3(il ] . Here we will apply it to a three-photon superposition 
state for the demonstration of three-photon de Broglie 
wave length without the NOON state. 

The NOON-state projection measurement scheme for 
three-photon case is sketched in inset (b) of Fig.l. 
In this scheme, the input field is first divided into 
three equal parts. Then each part passes through a 
phase retarder that introduces a relative phase differ- 
ence of 0, 27r/3, 47t/3 respectively between the H- and V- 
polarization. The phase shifted fields are then projected 
to 135° direction by polarizers before being detected by 
A, B, C detectors, respectively. It was shown that the 
three-photon coincidence rate is proportional to 

'' {NOON 3 \§ 3 ) 



P 3 <x 



18 



(28) 



where \NOON 3 ) = (\3 H , 0y) - |0 H , 3 V »/V2 and |$ 3 ) = 
co\3 H ,0v) +cx\2hAv) +c 2 \1 h ,2 v ) +c 3 |0tf,3y). Note 
that since |2, 1), |1, 2) are orthogonal to the NOON-state, 
their contributions to P3 are zero. Assuming that |co| = 
I C3 j = c and there is a relative phase of ip between H and 
V so that c — c,c 3 — ce l3v , we obtain from Ea. l|28[) 



P3 oc — (1 - cos3</?), 



(29) 
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FIG. 4: Experimental result for the NOON state projec- 
tion measurement. The data is least-square-fitted to P<io(l + 
V3 cos3</9 with V3 = 84% after background subtraction. 



which is exactly in the form of Eq.J2J) with N — 3, show- 
ing the three-photon de Broglie wave length. 



B. Experiment 

From Sect.IIB, we learned that a state of \2jj, ly) can 
be produced with two parametric down-conversion pro- 
cesses. This state will of course give no contribution to 
the NOON state projection since it is orthogonal to the 
NOON state. On the other hand, we can rotate the state 
by 45°. Then the state becomes [20I l2lj 



2 H ,h 



(30) 



which has the NOON state component with c = v/3/8. 

Experimentally, the three-photon state of \2h, ly) is 
prepared in the same way described in Sect.IIB and 
shown in Fig. 2. Different from Sect.IIB, the polariza- 
tions of the prepared state are rotated 45° by HWP1 to 
achieve the state in Ea. l3Ufl . The phase shifter (PS) then 
introduces a relative phase difference ip between the H 
and V polarizations and HWP2 is set at zero before the 
NOON state projection measurement is performed [In- 
set (b) of Fig. 2]. As before, a four-photon coincidence 
measurement among ABCD detectors is equivalent to a 
three-photon coincidence measurement by ABC detec- 
tors gated on the detection at D, which is required for 
the production of \2h, ly). 

Four-photon coincidence count among ABCD detec- 
tors is registered in 200 sec as a function of the phase 



J 



ip (PS). The data after subtraction of background con- 
tributions is plotted in Fig. 4. It clearly shows a sinu- 
soidal dependence on <p with a period of 27r/3. The 
solid curve is a chi-square fit to the function of P4 = 
P4o[l + V 3 cos3(v3 + tpo)} with P 40 = 103/ 200sec and 
V3 = 0.84. The \ 2 of the fit is 24.3, which is comparable 
to the number of data of 25 indicating a good statistical 
fit. 

The less-than-unit visibility is a result of temporal dis- 
tinguishability among the three photons produced from 
two crystals. It can only be accounted for with a multi- 
mode model of the state given in Sect.IIC. Let's now 
apply it to the current scheme. 



C. Multi-mode analysis 

The input state is same as Eg. (11411 . But the field op- 
erators are changed to 



E A (t)= [E+(t)-e^E_(t)]/V6 + 
E B (t) = [E+(t) - e l ^E^(t)e^^]/VQ- 
E c {t) = [E+(t) - e^£L(i)e i47r / 3 ]/\/6- 



(31) 



with 



£+(*) = [E H (t)+E v (t)]/V2, 
E_(t) = [E H (t)-E v (t)]/V2, 



(32) 



where we omit the vacuum input fields and ip is the phase 
shift introduced by PS in Fig. 2. The field operator for 
detector D is same as Ea. lfTH|l . 

As in Sect.IIC, the four-photon coincidence rate is re- 
lated to the correlation function in Ea. (|18fl and we can 
first evaluate E A (ti)E B (t 2 )Ec(t3)E D (ti). With the field 
operators in Ea. (|31|l . we obtained 



E A {h)E B {t 2 )Ec(t 3 )E D (t 4 ) 
= (HHVai + HVHa 2 



VHHa 3 ) /12V12, (33) 



with 



ai = 1 + e l3v + 2e' 1 ( 2v+27r / 3 ) + 2e l ^ +i7r ^ 3 '> 
a 2 = 1 + e i3 ^ + 2e l ^ +i */V + 2 e l ^+ 27r / 3 ) (34) 
a 3 = l + e i3 ^ + 2e l2 ^ + 2e^, 

where the notations are same as in Eci. (|19f) and we used 
the identity 1 + e i27r/3 + e l47r/3 = 0. As before, we 
also drop five terms that give zero result when they 
operate on \^> 2 ). Now we can calculate the quantity 
EA{ti)E B {t 2 )Ec(t 3 )E D (t 4 )\y 2 ), which has the form of 



= {[G(tl,t2,t3,t4)+G(*2,*l,*3,*4)]ai + [G(h,t3,t 2 ,U)+G(t3,t 1 ,t2,t i )]a 2 



[G(t 2 , h, h, U) + G(t 3 , t 2 , h, U)] a 3 } | vac), (35) 



r 



where G{t\,t2 ,^3,^4) is given in Ea. (|21[l . After the time 
integral, we obtain 



P 4 (NOON) 



72 



(l + V 3 cos3<p) (36) 



with 



v 3 (ivocw) 



A + 2£ 
2A + £ : 



(37) 



where .4 and £ are given in Eas. i|25l I26|) . Note that the 
terms such as cos 2<p, cos ip are absent in Ea. (|37(l even 
in the non-ideal case of £ < A. This is because of the 
symmetry among the three detectors A, B, C involved in 
the three-photon NOON state projection measurement. 
When spatial mode mismatch is considered, the visibility 
is changed to 



V 3 (NOON) = v{ 



3-4 



2£ 



2A + £' 



(38) 



With v\ and the quantity £/A obtained in Sect.IIC, we 
have V 3 (NOON) = 0.85, which is close to the observed 
value of 0.84 in Sect.IIIB. 



IV. SUMMARY AND DISCUSSION 

In summary, we demonstrate a three-photon de Broglie 
wavelength by using two different schemes of projection 
measurement without the need for a hard-to-produce 
NOON state. Quantum interference is responsible for 
the cancellation of the unwanted terms. The first scheme 
by asymmetric beam splitters targets specific terms while 
the second one by NOON state projection cancels all the 
unwanted terms at once. We use a multi-mode model to 
describe the non-ideal situation encountered in the exper- 
iment and find good agreements with the experimental 
results. 

Since the scheme by asymmetric beam splitters is only 
for some specific terms, it cannot be easily scaled up to 



arbitrary number of photons although the extension to 
the four-photon case is available. The extension of the 
scheme by NOON state projection to arbitrary number 
of photons is straightforward. In fact, demonstrations 
with four and six photons have been done with simpler 
states [l6lll7| . 

On the other hand, the scheme of NOON state pro- 
jection need to divide the input fields into N equal parts 
while the scheme with asymmetric beam splitters requires 
less partition. So the latter will have higher coincidence 
rate than the former. In fact, Fig. 3 and Fig. 4 show a 
ratio of 4 after pump intensity correction. This is consis- 
tent with the ratio of 4.8 from Eqs . ESI> when £ = A. 
The difference may come from the different collection ge- 
ometry in the layout. 

The dependence of the visibility in Eas. (|27l I37fl on the 

quantity £/A reflects the fact that the interference ef- 
fect depends on the temporal indistinguishabilit y o f the 
three photons. From previous studies |15l l23l l3ll 133 . l33| . 
we learned that the quantity £ /A is a measure of indis- 
tinguishability between two pairs of photons in paramet- 
ric down-conversion. In our generation of the |2#,ly) 
state, one of the H-photon is from another pair of down- 
converted photons. So to form an indistinguishable 
three-photon state, we need pair indistinguishability, i.e., 
£/A^l. 



Acknowledgments 

This work was funded by National Fundamental Re- 
search Program of China (2001CB309300), the Innova- 
tion funds from Chinese Academy of Sciences, and Na- 
tional Natural Science Foundation of China (Grant No. 
60121503 and No. 10404027)). ZYO is also supported 
by the US National Science Foundation under Grant No. 
0245421 and No. 0427647. 



[1] J. Jacobson, G. Bjork, I. Chuang and Y. Yamamoto, 

Phys. Rev. Lett. 74, 4835 (1995). 
[2] W. Heisenberg, Zeitschr. f. Physik 43, 172 (1927). 
[3] M. J. Holland and K. Burnett, Phys. Rev. Lett. 71, 1355 

(1993). 

[4] J. J. Bollinger, W. M. Itano, D. J. Wineland, and D. J. 

Heinzen, Phys. Rev. A54, R4649 (1996). 
[5] Z. Y. Ou, Phys. Re. Lett. 77, 2352 (1996); Phys. Rev. A 

55, 2598 (1997). 
[6] A. N. Boto, P. Kok, D. S. Abrams, S. L. Braunstein, 



C. P. Williams, and J. P. Dowling, Phys. Rev. Lett. 85, 

2733 (2000); P. Kok, A. N. Boto, D. S. Abrams, C. P. 

Williams, S. L. Braunstein, and J. P. Dowling, Phys. Rev. 

A63, 063407 (2001). 
[7] Z. Y. Ou, X. Y. Zou, L. J. Wang, and L. Mandel, Phys. 

Rev. A 42, 2957 (1990). 
[8] J. G. Rarity, P. R. Tapster, E. Jakeman, T. Larchuk, R. 

A. Campos, M. C. Teich, and B. E. A. Saleh, Phys. Rev. 

Lett. 65, 1348 (1990). 
[9] H. F. Hofmann, Phys. Rev. A 70, 023812 (2004). 



[10] F. Shafiei, P. Srinivasan, and Z. Y. Ou, Phys. Rev. A 70, 
043803 (2004). 

[11] B. Liu and Z. Y. Ou, Phys. Rev. A 74, 035802 (2006). 
[12] P. Walther, J. W. Pan, M. Aspelmeyer, R. Ursin, S. 

Gasparoni, and A. Zeilinger, Nature (London) 429, 158 

(2004). 

[13] M. W. Mitchell, J. S. Lundeen, and A. M. Steinberg, 

Nature (London) 429, 161 (2004). 
[14] H. B. Wang, and T. Kobayashi, Phys. Rev. A 71, 

021802(R) (2005). 
[15] F. W. Sun, Z. Y. Ou, and G. C. Guo, Phys. Rev. A 73, 

023808 (2006). 

[16] K. J. Resch, K. L. Pregnell, R. Prevedel, A. Gilchrist, 

G. J. Pryde, J. L. O'Brien, and A. G. White, 

|quant-ph/0511214| 
[17] F. W. Sun, B. H. Liu, Y. F. Huang, Z. Y. Ou, and G. C. 

Guo, Phys. Rev. A 74, 033812 (2006). 
[18] K. Sanaka, K. J. Resch, and A. Zeilinger, Phys. Rev. 

Lett. 96, 083601 (2006). 
[19] C.K. Hong, Z.Y. Ou, and L. Mandel, Phys. Rev. Lett. 

59, 2044 (1987). 
[20] Z. Y. Ou, C. K. Hong, and L. Mandel, Opt. Commun. 

63, 118 (1987). 

[21] R. A. Campos, B. E. A. Saleh, and M. C. Teich, Phys. 



S 



Rev. A 40, 1371 (1989). 
[22] L. Mandel and E. Wolf, Optical Coherence and Quantum 

Optics, (Cambridge University Press, New York, 1995). 
[23] B. H. Liu, F. W. S un, Y. X. Gong, Y F. Huang, Z. Y. 

Ou, and G. C. Guo, quant-ph/0606118 
[24] R. Ghosh, C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. 

Rev. A34, 3962 (1986). 
[25] Z. Y. Ou, Quan. Semiclass. Opt. 9 599 (1997). 
[26] W. P. Grice and I. A. Walmsley, Phys. Rev. A 56, 1627 

(1997). 

[27] D. Branning, W. P. Grice, R. Erdmann, and I. A. Walm- 
sley, Phys. Rev. Lett. 83, 955 (1999) 

[28] O. Kuzucu, M. Fiorentino, M. A. Albota, F. N. C. Wong, 
and F. X. Kartncr, Phys. Rev. Lett. 94, 083601 (2005). 

[29] Z. Y. Ou, quant-ph/0601118 

[30] G. Y. Xiang, Y. F. Huang, F. W. Sun, P. Zhang, Z. Y. 

Ou, and G. C. Guo, Phys. Rev. Lett. 97, 023604 (2006). 
[31] Z. Y. Ou, J. K. Rhee, and L. J. Wang, Phys. Rev. A 60, 

593 (1999). 

[32] K. Tsujino, H. F. Hofmann, S. Takeuchi, and K. Sasaki, 

Phys. Rev. Lett. 92, 153602 (2004). 
[33] Z. Y. Ou, Phys. Rev. A 72, 053814 (2005). 



